Vicinal surface growth: bunching and meandering instabilities 
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The morphology of a growing crystal surface is studied in the case of an unstable two-dimensional 
step flow. Competition between bunching and meandering of steps leads to a variety of patterns 
characterized by their respective instability growth rates. The roughness exponent is shown to go 
from 1/2 to 1, between the pure bunching to the meandering regimes. Using numerical simulations, 
we observe that generically, a transition between the two regimes occurs. We find surface shapes 
and roughness time evolution in quantitative agreement with experiments. 
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The interplay of instability and nonlinearity is central 
to the formation of spatiotemporal patterns in physical 
systems driven away from equilibrium [1]. The typical 
scenario for the appearance of a periodic structure in a 
fluid, a spiral in an excitable medium, or a specific front 
velocity in a reaction-diffusion system, start with the de- 
velopment of a linear instability, which fixes the char- 
acteristic scale, followed by some nonlinear saturation 
mechanism, which fixes the characteristic similarity ex- 
ponents. This scenario can drastically change in the pres- 
ence of two simultaneous linear instabilities, because of 
the competition between different time and length scales. 

The step flow [2, 3] during molecular beam epitaxy of 
a vicinal surface is subject to different kinds of insta- 
bilities such as bunching (^-direction) [4] and meander- 
ing (^-direction) |5] of steps. The intrinsic anisotropy of 
the vicinal surface leads to different diffusion and attach- 
ment mechanisms and ultimately, to coarsening power 
laws with different scalings according to the directions 
along, or perpendicular to the steps [6, 7J. Under the 
simplest conditions of homoepitaxy, bunching and mean- 
dering develop together, as for instance in copper |8l4l0j, 
and (001) silicon surfaces [Ill4l5l|, and are at the origin 
of a rich variety of surface patterns depending on the 
growing conditions. 

In this Brief Report, we propose a "vicinal equation" to 
describe, in a long wavelength and small amplitude ap- 
proximation [1611, th e simultaneous development of both 
instabilities [14|, [l5| . The evolution of the surface height 
z = Ft — mx + h(x, y, i), with F the deposition flux, m 
the mean slope and h the surface shape, is given by, 

ht=Lh- V 2 |V/i| 2 - AV • (h x Vh) , (1) 

where the operator L is defined by the expression 

L={a + \k x -kl-k 2 y )kl + {b-k 2 y )k 2 y , 

in Fourier space. The system's behavior is determined 
by the intrinsic anisotropy of the vicinal surface, and de- 
pends on the respective characteristic growth times, fixed 



by two parameters, a and 6, proportional to the bunching 
and meandering amplification rates, respectively. The 
first term, in L, contains the linear effects describing the 
instabilities, step flow dispersion and diffusion; the sec- 
ond term is the coarsening nonlinearity and the last one 
the stepping nonlinearity. 

Equation (pQ) amounts to a generalization of the 
anisotropic evolution of a vicinal surface, unstable with 
respect to meandering but stable to bunching, of Ref. jgj]. 
The essential difference is the presence of a new nonlin- 
earity, the last term in (pQ), proportional to the coupling 
parameter A, that breaks the symmetry x O — x, and 
allows us to take into account the mechanisms leading to 
step bunching. It can be viewed as a simple small gra- 
dient expansion of step flow driven by a non-equilibrium 
current [171 ], 
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where S = 1+ — L is the difference between the up step 
/ + and down step L attachment lengths; in this case one 
obtains A = F5/m 3 , and bunches form if S < 0. 

We may consider Eq. (PQ) as a phenomenological model 
of the vicinal surface growth; the actual form of the three 
independent parameters (a, 6, A) depends on the physi- 
cal details of the meandering and bunching instabilities. 
Its special form assumes that the essential mechanisms 
of these instabilities are kinetic, as appropriated for ho- 
moepitaxy conditions. Throughout this paper we will 
use units based on the equilibrium terrace length Zq, and 
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FIG. 1. (Color online) Surface patterns viewed from bottom 
—h(x,t); bunches, meanders and cellular patterns, from left 
to right. 
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FIG. 2. (Color online) Bunching instability (left) a = 1.0, b = 0.1, and meandering instability (right) a = 0.1, b = 1.0, in the 
weak stepping nonlinearity regime (A = 0.2). Roughness w(t) (top) and length scales \ X) X y (bottom). Snapshot of the surface 
shape at t = 400 (inset). The similarity exponents (3 and a are shown to schematically indicate the different regimes (insets in 
log- log scale). 



time Iq/DoCo, with adatom concentration Co and diffu- 
sion coefficient Do, conveniently rescaled to throw out 
non relevant dependent parameters [6]. Therefore, com- 
parison with experiments can be done by choosing length 
and time scales to fit a given reference state. 

For different sets of parameters a, b and A, meanders, 
step bunches or cellular patterns may be obtained, as 
illustrated in Fig. [U The snapshots show the surface 
shape, obtained after the evolution of an initially homo- 
geneous state perturbed by small amplitude noise, viewed 
from the bottom, to better visualize the strong gradients 
formed in both directions. 

We can distinguish two limiting one-dimensional ap- 
proximations according to the relevant unstable direc- 
tion. In the bunching dominated regime the effective 
equation is, 

h t ~ -ah xx - \{h 2 x ) x (3) 

(we use subscripts to denote partial derivatives). Putting 

h(x,t) = t^H{x/t° ix ), (4) 

one obtains a x = /3 = 1/2, where the exponents a x and /3 
govern the characteristic bunching length and the rough- 
ness time scales, respectively. In the meandering regime 



we have, 

kt ~ -bkyy - (h 2 y )yy (5) 

and 

h(y,t)=t^H(y/t a y) (6) 

which gives /? = 1 for the roughness exponent and 
Oiy = 1/2, for the meander length scale. In addition, 
the stable direction is governed by the bi-quadratic equa- 
tion h t ~ — V 4 /i, leading to a characteristic spreading 
exponent 1/4, in the x-direction for bunching, and in the 
^-direction for meandering. Therefore, one may predict 
that the roughness exponent j3 should go from 1/2 in 
the pure bunching, to 1 in the meandering regime. We 
performed a series of numerical computations to obtain 
the asymptotic behavior of the surface roughness, and to 
investigate the pattern formation in the presence of the 
anisotropic instabilities. Our goal is to compare these re- 
sults with experiments of Si (001) homoepitaxy pj]-[T3j. 

To illustrate the two limiting regimes, we choose two 
sets of linear parameters: (i) a = 1.0, b = 0.1 for the 
bunching case, and (ii) a = 0.1, b = 1.0 for the meander- 
ing case. The actual behavior of the system also depends 
on A, we distinguish between a weak nonlinearity case 
A = 0.2, and a strong one A = 1.0. The system size is a 
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FIG. 3. (Color online) Bunching instability (left) a = 1.0, b — 0.1, and meandering instability (right) a = 0.1, b = 1.0, in 
the strong stepping nonlinearity regime (A = 1.0). Roughness w(t) (top) and length scales A^, X y (bottom). Snapshot of the 
surface shape at t — 160 (inset). The similarity exponents /3 and a are shown to schematically indicate the different regimes 
(insets in log- log scale). 



256 x 256 square (we use a pseudo-spectral algorithm with 
512 2 collocation points). Results in these two regimes are 
displayed in Figs. [2] and [3] We plot the roughness, 



W(t) = y/(h(t)*) - (h{t)Y 



(7) 



where the average is over the whole surface, and the char- 
acteristic lengths, 



(A x 2 , A 2 ) 



Jdk(klP y ) \h k f 
J dk \h k \ 2 



(8) 



where h k = h k (t) is the two-dimensional Fourier trans- 
form, together with the corresponding snapshots of the 
surface shape. The lengths X y (t) and X x (t) are the char- 
acteristic scales in the meander and the step gradient 
directions, respectively. 

In accordance with (HJ) and (|6]), we may expect the 
asymptotic behaviors: w ~ t 13 and X x ~ t% X y ~ t a y . 
One may think in particular that, in the bunching case, 
the roughness evolution actually observed would start, 
after a transitory dependence on the initial state, with 
a w ~ t x l 2 growth (dominant at small times), followed 
by a w w t at larger times. This is effectively the case 
when the stepping nonlinearity is dominated by coarsen- 
ing (Fig. [2l top- left). Indeed, after an initial evolution 



characterized by the bunching instability with j3 — 1/2, 
the coarsening regime is established at long- times, lead- 
ing to a roughness evolution with /3 = 1. 

However, at long times, depending on the value of A, a 
new regime may arise. It depends on the relative strength 
of the coarsening nonlinearity and the stepping nonlin- 
earity. These become comparable at the scale: 



V 2 |V/i| 2 - AV • (h x Vft) . 



Assuming A^ 



X T 



I <C L to ensure that two- 



dimensional effects are important, this can be satisfied 
if 



1 
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where L is the system size. In the step bunching unstable 
flow, when Xl(t) becomes greater than one, the stepping 
nonlinearity overcomes the coarsening one. In this case 
the roughness would be characterized by an exponent 
f3 = 1/2, as found in the strong bunching regime (Fig. [3j 
top- left). We then observe an initial regime dominated by 
coarsening (with /3 = 1), followed by a regime dominated 
by the stepping nonlinearity (with /3 = 1/2). It is worth 
noting that the change in regime takes place around t w 
300, when the X x scale overcomes the A^ scale; in the 




weak stepping nonlinearity regime of Fig. [2l this never 
happens. 

In the case where b > a the meandering instability 
dominates (Figs. [2] and [3j right). For both cases, weak 
and strong stepping nonlinearity, the X x (t) characteris- 
tic length rapidly increases, signaling that a quasi one- 
dimensional state sets in. The roughness follows a power 
law with j3 = 1. The phenomenology of the meander 
coarse graining is similar to the one described by the weak 
nonlinear amplitude equation for the Bales- Zangwill in- 
stability d, [Hi. It is interesting to note the meander 
morphology is modified by the presence of a bunching 
instability, although weak, that accelerates the homoge- 
nization along the steps, eliminating the gradients in this 
direction. 

In the bunching regime a > 6, for weak (but unsta- 
ble) meandering and step nonlinearity, the shape of the 
surface consists on a succession of bunches with some 
dislocation type defects. It compares satisfactorily with 
the atomic force microscopy (AFM) measurements of a 
silicon vicinal (001) surface [13|, as shown in Fig. HJ The 
observed quasi-periodic pattern of step bunches and the 
time scales are well reproduced by the numerical model 
(PQ). In particular, the observed experimental evolution 



of the roughness, that the model predicts to start with a 
first asymptotic with /3 = 1/2 afterward increasing up to 
/3 = 1, actually shows a breaking between a slow initial 
grow that later becomes steeper; a fit over the full time, 
mixing the two regimes, gives an intermediate exponent 
Pexp — 0.53. 

In this Brief Report we presented a weak nonlinear 
model of the evolution of a vicinal surface subject si- 
multaneously to bunching and meandering of steps. We 
found that depending on the relevant nonlinear mech- 
anism, coarsening or stepping, different long time evo- 
lutions of the surface roughness are possible. We also 
showed that the presence of a weak, subdominant, bunch- 
ing instability, modifies the morphology of the meanders, 
leading to a quasi one-dimensional pattern; the same is 
true for the bunching dominant regime, when coarsening 
is the relevant nonlinear effect. In conclusion, the phe- 
nomenological vicinal surface equation is able to repro- 
duce the observed morphology and roughness evolution 
of homoepitaxial growth, and provides an useful frame- 
work for investigating patterns in anisotropic unstable 
systems. 
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